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ABSTRACT 

Serendipitous stellar occultation search is so far the only way to detect the existence 
of very small, very dim, remote objects in the solar system. To date, however, there 
are only very few reported detections for trans-Neptunian objects (TNOs) in optical 
bands. In the X-ray band, with the RXTE/PCA data of Sco X-l taken from June 
2007 to October 2009, we found one possible X-ray occultation event. We discuss the 
veracity and properties of this event, and suggest upper limits to the size distribution 
of TNOs at hectometer size and of Main-Belt Asteroids (MBAs) at decameter size. 

Key words: occultations - Kuiper Belt - Solar system: formation - stars: neutron 
- X-rays: binaries. 



1 INTRODUCTION 

Population properties of small solar-system bodies, includ- 
ing main-belt asteroids (MBAs) and trans-Neptunian ob- 
jects (TNOs), among others, carry information of the early 
solar system when planets were being formed and of the 
physics in the dynamical and collisional histor y of our solar 
system (see, for example. iKenyon et al, iBottke et al.l 

(20051 ): io'Brien fe Greenbergj (|2005l ): Ichend (120041 1 ). Our 
knowledge of them, however, is far from being complete, 
particularly for those very small ones. Direct observation is 
so far achieved only for TNOs larger than deca-kilometers 
and for MBAs larger than about 400 meters. To pin 
down their size distribution at the small size end, search- 
ing for serendipitous stell ar occultation events was pro- 
posed to be a possible way (iBailev 19761: Brown fc Webster] 



Il997i : iRoaues fc Moncuquetoud : Icoorav fc Farmer! l2003h . 
Such se arch has been being conducted mainly in optical 
band s (jBianco et al.l 201ol : iBickerton. Kavelaars fc Welch! 
120081: Rogues et al.l 12003 ) with few r eport ed detections 
( Rogues et alj l200d : ISchlichting et all 120091 1. In the X- 



ray band, in which occultation by even smaller objects 
may be more easily detected because of a smaller Fres- 
nel scale, putative occultation events, allegedly caused by 
100-meter size TNOs, were found in the 1996-2002 X- 
ray data of Sco X-l taken by the instrument Propor- 
tional Counter Array ( PCA) on boa r d Ro s si X-r ay Tim- 
ing Explorer (RXTE) (|Chang et all l200d . l2007h . Those 
events, however, were later found to be most likely re- 
lated to some instrumental dead-time effect of unknown na- 
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ture dChang et all 120071; I Jones et al.ll200Sl ; [Liu et all W0i : 
iBlocker. Protopapas fc Alcockl 120091 ). To get rid of instru- 
mental effect contamination, new RXTE/PCA observations 
of Sco X-l with a newly designed data mode to record de- 
tailed information of each PCA-detected high energy event 
have been being conducted since June 2007. It is hoped that 
those high energy events, which are called 'very large events' 
(VLEs) in RXTE/PCA terminology, can provide clear sig- 
nature of possible instrumental effects. 

In this paper we report the result we obtained from 
analyzing the new data taken from June 2007 to October 
2009, which in total is about 240 ks. We found that VLEs, 
classified into different types according to the number of 
triggered anodes recorded in the data, can be used as an 
indicator of instrumental effects. In the 240-ks data, only 
1 dip event is found to be unrelated to VLEs. We discuss 
the property of this dip event with occultation diffraction 
pattern fitting to its light curve. We then estimate upper 
limits to the size distribution of TNOs at hectometer size 
and of MBAs at decameter size, which are unattainable by 
any other means known to date. 



2 DIP EVENTS AND THE RXTE/PCA VLES 

The new RXTE observation of Sco X-l started from June 
13, 2007, and is still going on (RXTE Obsld 93067). Tech- 
nical details of RXTE instr umentation can be f ound in the 
RXTE web site and also in ljahoda et al.l ([2006). To search 
for millisecond dip events, we a dopted the same proce- 
dure employed in previous studies ijChang et al . 2006, 2007; 
lLiu et al.ll2008l h We examined the light curves of Sco X-l 
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Type 


number of triggered anodes 


average count rate 


A 


no 


1.73 ±0.29 


B 


all 


34.5 ± 11.2 


C 


more than one but not all 


44.8 ±9.80 


D 


only one 


9.44 ±2.27 



Table 1. VLE types and their average count rates (in units of 
counts per second per PCU). The number of PCUs which are 
on during the observation varies from time to time. The total 
length of the data employed in this analysis is 240 ksec, and is 
970 ksec-PCU when the number of PCU on is taken into account. 



binned in different bin sizes from 1 ms to 10 ms. The de- 
viation in photon counts of each bin in the light curve was 
determined by comparing its count number with the average 
counts and the count variance in an 8-s running window. In 
the 240-ks data that we use for this analysis, 39 'significant' 
and 253 'less significant' dips are found. To compare with 
our earlier works, we define 'significant' as those dips found 
in our searching algorithm with negative deviation larger 
than 6.5 a and 'less significant' as that between 5 a and 6.5 
a. The choice of (minus) 6.5 a in previous works was to set 
a random probability at 0.001 for the a nalys is of the original 
320-ks data employed in I Chang et alj ((2006) , and (minus) 5 
a is about the level below which the deviation distribution 
obviously shows excess compared to a normal distribution. 
We keep the same choice of the deviation level as in earlier 
works to have a better c o mparison. Read ers are referred to 
IChang et all l|2006l . l2007l ); lLiu et al l (120081 ). for details of our 
searching algorithm. To investigate the connection between 
these dip events and possible instrumental effects, we exam- 
ined detailed information of the recorded VLEs as described 
below. 

RXTE/PCA VLEs are those events which deposit more 
than 100 keV in any one of the 6 active xenon anodes and 
the propane anode of one of the five identical proportional 
counter units (PCUs) of PCA. A VLE saturates the pre- 
amplifier and causes ringing in the signal chain when the 
amplified signal is restored to the baseline. Longer system 
dead time is set for VLEs, usually at 50 fis for observations of 
Sco X-l. The actual duration in which the signal chains are 
affected depends on the pulse height of the saturating event. 
In default setting of RXTE/PCA observations, VLE count 
rates are recorded in a standard data mode with 125- ms time 
resolution. In the new observation of Sco X-l (RXTE Ob- 
sID: 93067), to explore possible instrumental effects which 
may cause the millisecond dips in the PCA light curve of Sco 
X-l, a particular data mode was designed to retain detailed 
information of each individual VLE. In this data mode, the 
event epoch is recorded with 125- /^s resolution and the iden- 
tification of anodes triggered by the VLE is also recorded. 
With this new data mode, we found that some VLEs are 
very unusual in the way that no anode is recorded as hav- 
ing been triggered. We therefore classify VLEs into 4 types 
according to the number of triggered anodes as listed in Ta- 
ble [I] in which the averaged count rates (in units of counts 
per second per PCU) of these different VLEs are also listed. 
The count rate for all VLEs is 90.5T19.2 counts per second 
per PCU. It is clear that Type-A VLEs are rare. 

To study the association of the dip events that we found 



Group associated VLE type significant less significant 







dips 


dips 


A 


Type A 


34 


125 


B 


Type B, no Type A 


2 


94 


C 


Type C, no Type A, B 


2 


23 


D 


Type D only 





3 


E 


no VLEs 


1 


8 



Table 2. Number of dips in different groups. 



with VLEs, we classify the dip events into 5 groups according 
to the type of VLEs that were detected within the time 
interval defined by the dip duration and one duration before 
the dip epoch. The number of dips in different groups is 
shown in Table [2] It is obvious that dip events are strongly 
related to Type-A VLEs. Among the 39 significant dips 34 
are associated with Type-A VLEs, which are rare among 
all VLEs. It is not clear yet why a VLE is recorded in the 
data with no anode triggered. VLEs are usually caused by 
high energy particles. It may be possible that the incident 
particle is very energetic and an ionization avalanche or a 
secondary-particle shower is created so t hat the instrumen t 
suffers from a long sequence of dead time (|jones et al ]|2008l ). 
For some unknown reasons all the VLE triggering flags are 
turned off. Although not yet fully understood, we conclude 
that Group-A dip events are instrumental. 

To investigate whether dip events in other groups are 
also instrumental, we examine the consistency between the 
number of less significant dips and that due to random fluc- 
tuation. The deviation distribution of all the time bins in 
the light curve should be of a Poisson nature modified by 
the dead-time and coincidence-event effects if all the fluctu- 
ations are random and there is no occultation event. One ex- 
ample is shown in Figure Q] for a search with 2- ms time bins. 
Similar figur e s for earlie r data can be found in lChang et al.1 
(|2006l . 120071 '); iLiu et all (|2008T ). From Figure [T] we can see 
that, although Gaussian distribution is a good approxima- 
tion, the deviation distribution of the data is about a factor 
of 2 ~ 3 lower than the Gaussian one at — 5cr. In our cur- 
rent discussion we mainly consider the search with 2-ms bins 
since most of the less significant dips have a duration of 2 
ms. If a Gaussian distribution is assumed, one expects to 
have 34.4 bins with photon counts less than 5 a below the 
average for 2-ms bins in the 240-ks data. Since some VLEs 
happen at the same time (within the 125ps time resolution) 
or within the same 2-ms bin, to estimate the distribution 
of these 34.4 bins in different dip groups, we first count the 
number of bins without any VLEs. The total amount of 
those bins is about 64.5 ks in the 120-ks data. Therefore, 
18.5 of the 34.4 bins should be associated with Group E. 
For the other dip groups, we distribute the remaining 15.9 
bins according to their related VLE count rates. The result 
is that, assuming a Gaussian distribution, there should be 
0.3, 6.1, 7.9, 1.7, and 18.5 bins associated with Group A 
to E respectively. We note that this is in fact not strictly 
correct, because as mentioned above some VLEs happen in 
the same 2-ms time bin and we assign dip events into differ- 
ent groups with highest priority for Type-A VLEs, that is, 
there can be Type-B VLEs associated with Group-A dips. 
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Figure 1. Deviation distribution of the Sco X-l RXTE/PCA 
light curve based on the 240-ks data taken from June 2007 to 
October 2009 that we employed in this work. Thick histograms 
are from the data with 2-ms time bins, and thin histograms are the 
corresponding Gaussion distribution with the same total number 
of bins. The upper panel is for all the time bins, and the lower 
panel is for only those bins within which there are no VLEs. In 
the lower panel, most of those bins with large negative deviation, 
except for the one at about —7.1c, arc related to dip events of 
Group A-D, although within the 2-ms bin itself there is no VLE. 



Nontheless, the information here is enough for us to draw 
conclusions. 

Comparing the numbers of less significant dips in differ- 
ent groups with those derived from a Gaussian distribution, 
and further taking in to account the reduction by a factor of 
2^3, one can see that the numbers of less significant dips of 
Group A, B, C and D are not consistent with random fluc- 
tuation. There is a tendancy that the excess of dip numbers 
decreases from Group A to Group D. Although details are 
not yet fully understood, it seems clear that Type-A VLEs 
are the most energetic and Type-D VLEs are the weakest. 
Because of the association with VLEs, all the Group A-D 
dips that are not due to random fluctuation are most likely 
instrumental. On the other hand, the number of less sig- 
nificant dips of Group E is roughly consistent with random 
fluctuation. We have no indication of instrumental effects for 
Group-E dips, which are not associated with any VLEs. The 
significant dip in Group E, which is found in the search with 
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Figure 2. Light curves of Event El. The upper panel shows the 
light curve 100-ms before and after Event El in 0.5-ms bins. The 
lower panel is a zoom-in view of Event El. The thick histogram 
is the observed light curves in 0.5-ms bins, the thin curve is the 
computed diffraction pattern with best-fit parameters, and the 
thin histogram is the binned light curve of the computed one. The 
level of one-er uncertainty in the observed light curve is indicated 
with a vertical bar. The average number of counts in a 0.5-ms bin 
is 13.1. 

2-ms bins, has a deviation of —7.1a, corresponding to a ran- 
dom probability of less than 7.5 x 10 -5 in the 240-ks data. 
We denote this dip event as Event El, whose light curve 
is plotted in Figure [21 together with the best-fit diffraction 
pattern described in the next section. 



3 DISTANCE AND SIZE OF THE OCCULTING 
BODIES 

Event El, a non-instrumental millisecond dip event, is most 
likel y an occultation event cau sed by objects in the solar sys- 
tem l|Chang et al JI2006L l2007h ■ The observed light curve, in- 
deed, can be well described by shadows of occultation when 
diffraction is properly taken into account. In the following 
we explain how we conduct the diffraction pattern fitting 
and how we estimate the distance and size of the occulting 
body with fitting results and assumed orbital parameters. 

3.1 Fitting the occultation light curve with 
diffraction patterns 

In our study we assume the occulting body is a spherical one 
and the background star, Sco X-l, is a distant point source. 
We also take into account the actual RXTE/PCA observed 
Sco X-l spectrum, which peaks at 4 keV. We note that the 
shadow diffraction patterns are all the same for a given ratio 
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of the occulting body diameter to the Fresnel scale when the 
length is expressed in units of the Fresnel scale. Hereafter 
we denote that ratio with a. To describe an occultation light 
curve, 4 more parameters are needed. One is the impact pa- 
rameter ft, which is defined as the closest distance of the 
crossing path to the shadow center in units of the occulting 
body radius. Another one is the transverse relative speed ut 
of the observer with respect to the shadow. The third is the 
epoch of the central crossing, to, and the fourth is the noma- 
lization TV. For each occultation light curve, to find a best 
fit in the 5-parameter space (a, P, u t , to, TV) is not practical. 
We therefore search for the 'best' fit in the following way. 
We fix P at certain different values, e.g., from to 1 with 
an increment step of 0.1. For a fixed j3, we compute theo- 
retical occultation light curves over a range of a. For each 
computed light curve of a given a, u t over a certain range is 
applied to convert the length of the computed pattern into 
time and the computed curve is folded into a binned curve 
to compare with data. 

In this study, we use data curves with 0.5-ms bins. The 
normalization TV of the computed curve is set to be the av- 
erage counts per bin of the data curve in the fitting window, 
which we choose to be 12.5 ms, that is, 25 bins. With given 
P, a, u t , and TV, we adjust to with a step size of 0.05 ms, 
which is 0.1 time bin, to look for the smallest \ 2 between 
the binned curve and data curve. We then assign this \ 2 
to be the fitting result of this set of parameters /3, a, and 
u t . The 'best' fit is determined by the minimum \ 2 m the 
parameter space of P, a, and u t as described above. 

The best fit result for Event El is that x 2 = H-7 
(with 20 degrees of freedom), P = 0.0, a = 5.2_2.6 and 
ut = 2700_i45o Fresnel scale per second. The uncertainty 
quoted above is at the 3-er level in the a-u t space for a fixed 
P. Here we report only the lower uncertainty level, since a 
combination of large u t and large a gives acceptable fits; see 
Figure[3] The reason is that the light curve of Event El does 
not show clear side lobes. Diffraction patterns of large ut and 
a therefore cannot be discriminated. Those patterns corre- 
spond to occultation caused by very nearby objects (see the 
next subsection). The uncertainty in a quoted here is the 
value corresponding to the uncertainty in u t , rather than 
the a extreme values on the 3-cr level contour. It is because 
that our interest in a is for estimating the diameter of the 
occulting body, which depends on both u t and a in a non- 
trivial way, and the distance estimate depends on ut only. 
The P dependence of the diffraction pattern fitting is weak. 
The best-fit x 2 value (11.7 for 20 d.o.f.) is very low. It only 
shows that the light curve can be described by a diffraction 
pattern, but at the same time provides less power in dis- 
criminating fitting parameters. This low \ 2 value are due 
the following two factors. One is the low count number in 
each bin, which gives a relatively large error bar, and the 
other is the 0.5-ms bin size, which is large so that details 
of a diffraction pattern are smeared out. These two factors 
make it easy to obtain a low \ 2 value fitting. 

3.2 Distance and size determination assuming a 
circular orbit 

Based on the best estimated value of the relative transverse 
speed between RXTE and the shadow obtained from diffrac- 
tion pattern fitting and the knowledge of RXTE's velocity 
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Figure 3. \ 2 contours of the diffraction pattern fitting in the 
a-ut parameter space for Event El. The cross symbol marks the 
location where the smallest \ 2 value is found. 



at the event epoch (given in Table [3]), one may infer the dis- 
tance to the occulting body by assuming its orbital parame- 
ters. This assumption is needed because of the size-distance 
degeneracy, as one can see that Ut can be determined from 
diffraction pattern fitting only in units of Fresnel scale per 
second and the Fresnel scale is distance dependent. Since 
most of MBAs and TNOs have low-eccentricity orbits, for 
simplicity and as an approximation, we consider only circu- 
lar orbits in this study, but allow different orbital inclination. 
In Figure [4] the relative speed transverse to the Sco X-l di- 
rection between RXTE and the occulting body as a function 
of heliocentric distance of the occulting body is plotted, and 
so is the relative speed determined from diffraction pattern 
fitting, which is translated from the form of Fresnel scale 
per second into real speed as a function of heliocentric dis- 
tance of the occulting body. The intersection of these two 
curves gives the estimate of the heliocentric distance of the 
occulting body. From that distance, one may infer the di- 
ameter of the occulting body from the best estimated value 
of the diameter to Fresnel scale ratio a, also obtained from 
the diffraction pattern fitting. 

From Figure 21 considering the curve for the minimum 
inclination (the thick solid line in Figure [4]), we obtain that 
the heliocentric distance and diameter of the occulting body 
are r = 3.4±2 4 4 AU and D = 38±l? m. The uncertainty 
quoted here is based on the uncertainty of the diffraction 
pattern fitting. Because the fitting does not provide con- 
straint for high tit, we do not have corresponding constraint 
at the small distance end in Figure [4] In fact, a very nearby 
object of meter size and with speed at a few kilometers per 
second could produce a light curve like Event El. This cor- 
responds to high u t and high a because of the very small 
Fresnel scale at a very nearby distance. We therefore ar- 
tificially assign the lower uncertainty level to r and D in 
such a way to symbolically indicate such a situation. On the 
other hand, the curve for the largest relative speed (the thick 
dashed line) in Figure [4] at about 40 AU corresponds to an 
inclination angle of 172° (retrograde orbits). If we consider 
high-inclination orbits, within the 3-cr fitting uncertainty, 
the occulting body could be at about 40 AU (and therefor 
about 150 m in size). Of course, such orbits are rare. 
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Epoch 


RXTE location (km) 


RXTE velocity (km/s) 


MJD 54622.397997210 
(2008-06-05) 


X = -3.947533 x 10 7 
Y = -1.338443 x 10 s 
Z = -5.803026 x 10 7 


V x = +3.482881 x 10 1 
Vy = -1.043335 x 10 1 
V z = -9.082142 x 10" 1 



Table 3. Spacecraft information for Event El. The RXTE location and velocity at the event epoch are expressed in the celestial 
equatorial coordinate system. The X-axis is in the direction of the vernal equinox, the Z-axis is in the north celestial pole direction, and 
the solar system barycenter is at the origin. These values are obtained from RXTE housekeeping data and the ephemeris of the Earth 
(JPL Horizons On-Line Ephemeris System, http://ssd.jpl.nasa.gov/horizons.html). The direction of Sco X-l is R.A. = 16 h 19 m 55.07 s , 
Dec.= — 15°38'25.0" and its solar elongation at the event epoch is 169.3°. 



Because of the size-distance degeneracy and the weak 
constraints that diffraction pattern fitting to the Event-El 
light curve can provide, we are not able to draw firm conclu- 
sions on the distance and size of the occulting body causing 
Event El. The size-distance degeneracy may be resolved by 
the knowledge of the occulting-body orbit, which is usually 
not available and needs to be assumed, or by an array of de- 
tectors, which, although not yet available now, may provide 
information of the shadow size directly. In the following sec- 
tions we discuss what constraints we may have on the MBA 
and TNO size distributions based on the search over the 
240-ks RXTE/PCA data of Sco X-l. The result is plotted 
in Figure^ We note that, limited by the RXTE/PCA count 
rate of Sco X-l, which is about 2 x 10 4 counts per second per 
PCU, our dip-event search algorithm can detect occultation 
events caused by o bjects of size only down to one or two 
times Fresnel scale (|Chang et al.| [2007). 



4 THE MBA SIZE DISTRIBUTION 

Direct observation of the main belt asteroids has only been 
achieved down to the size of about 400 meters. Earlier sur- 
veys revealed a wavy size distribution of MBAs with a bump 
around 100 k m and a secon d one around 5 km (se e Fig. 4 in 
IChend (|2004n and Fig. 14 in lO'Brien fc Greenbersl (120051) for 
a compiled distributi on; see also iTedesco fc Desert! (|2002l ) 
for the ISO result, jlvezic et alj i|200ll ) for SPSS, and 
lYoshida fc Makamural i|2007l ) for Subaru survey). It appears 
that most of the results are consistent with a power index of 
-3 (5 - 40 km) and of -1.3 (0.4 - 5 km) for the cumulative dis- 
tribution. An estimate of the cumulative number of MBAs, 
based on the ISO result, is (1.2 ± 0.5) x 10 6 for MBAs with 
diameter larger than 1 km. 

To estimate the occultation event rate, we formulate the 
MBA size distribution as the following: 



N- 



>D 



(1) 



with £>! = 5 km, b 1 = 3, JVl = 1.5 x 10 J for 5 km < D < 
40 km, £> 2 = 1 km, b 2 = 1.3, N 2 = 1.2 x 10 6 for 0.4 km 
< D < 5 km, and D 3 = 0.4 km, b 3 = 2.5, iV 3 = 3.9 x 10 6 
for D < 0.4 km. In this formulation, the normalization is 
set according to the ISO result, and the power index 2.5 for 
the state of collisional equilibrium with a constant tensile 
strength is assumed for D < 0.4 km. The expected event 
rate can be estimated as 



He — - — 

t Ha 



(2) 



(AU) 



Figure 4. Estimate of the occulting body distance for Event El. 
The thick curves are the relative speed transverse to the Sco X-l 
direction between RXTE and the occulting body, whose orbit is 
assumed to be a circular one. The thick solid curve is for the case 
of a minimum orbital inclination. The thick dashed curves give the 
largest and smallest relative transverse speed when all the possible 
orbital inclination is considered. The lower thick dashed curve 
almost coincident with the thick solid one. In this computation, 
RXTE velocity in the solar system at the event epoch and the Sco 
X-l direction, which is 5.5 degrees north of the ecliptic, are both 
taken into account. The thin curves are the relative transverse 
speed determined from diffraction pattern fitting. The thin solid 
one is the best fit and the thin dashed one give the uncertainty 
at the 3-<r level. For Event El, diffraction pattern fitting does not 
provide good constraint to the higher end of the relative speed, 
because of lacking of side lobes. 

where Q T is the sky area swept by all MBAs in the time 
interval of r, and J1a is the sky area over which MBAs are 
distributed. We then have, for MBAs with diameter larger 
than 10 m, 



10m 
V 



A 2 AD 



= (^) r x 6.6 x 10 s km 
A 2 



(3) 
(4) 



where typical values of v, the relative speed, and A, the dis- 
tance from the earth to the MBA, need to be assumed, and 
about 99.3% contribution comes from those with diameter 
smaller than 0.4 km. Adopting A = 2 AU and v — 12.6 
km/s, we have 



Q T = 3.7 x 10" r s" 



(5) 



The ecliptic latitude distribution of about 13,000 MBAs 
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with r* < 21.5 is shown in Fig. 12 of llvezic et al.l l|200ll) . 
Considering that Sco X-l is at 5.5° north of the ecliptic, 
we adopt an equivalent sky coverage area as being spanning 
over 20° around the ecliptic, that is, 



• 360 x 20 x 



— V 

180/ 



2.2 



(6) 



The expected event rate, R a is therefore about 1.7 x 10 
s" 1 . 

If we consider that Event El is caused by an MBA, 
the measured event rate with one detection (Event El) is 
R m = (1/240 ks)« 4.2 x 10" 6 s" 1 . R m is about 250 times 
larger than _R C . However, if that is the case, the MBA causing 
Event El is about 40 m in size. The chance is low that one 
40-m object is found but no smaller ones, unless the number 
of MB As drops significantly for size smaller than 40 m. If so 
and still assuming a single power law for size smaller than 
400 m, to have one detection in 240 ks for MBAs down to 40 
m in size, the power index 63 will need to be 5.8. It is a very 
steep power law, as shown in Figure [5] and also requires the 
paucity of smaller MBAs. We consider this possibility not 
plausible. 

If Event El is not caused by an MBA, we have no detec- 
tion of MBAs in the 240-ks data. Our dip-event search algo- 
rithm can detect occultation events caused by objects of size 
only down to one or two times Fresnel scale (|Chang et al.l 
l2007ft . For the MBA distance, it is about 10 m. Based on the 
non-detection, we estimate the upper limit to the size distri- 
bution at the level of assuming one detection. As calculated 
above, this level will give a detection rate 250 times larger 
than R e . To have such a detection rate, if we still adopt a 
single power law for the size range between 10 m and 400 
m, the power index 63 needs to be 4.0, instead of 2.5. This is 
quite steep, although a power index larger than 2.5 is usu- 
ally e xpected. In the collisional e quibrium, one expects to 
have i|0'Brien fc Greenberd [20031 ) 



7 + s/3 
2 + s/3 



(7) 



where s is the power index of the strength-size relation. For 
&3 = 4.0, we have s = —2.25, which is larger than all the 
results reported in the literature to date. 



5 THE TNO SIZE DISTRIBUTION 

We describe the TNO size distribution in terms of the differ- 
ential surface density at ecliptic. More specifically, we plot 

versus s, where dn is the 



in this paper the quantity 



dn 
d log s 



differential number of TNOs per square degree of size 
between s and s + ds and s is the diameter. 



5.1 RXTE observations 

For a background point source, the event rate is 



N 
~T 



f s 2 I dn\ 
J S1 \ds) 



sv ds 



d? 



,180, 
x (— ) 



(8) 



where iV is the number of detected events (assuming a 100% 
detection efficiency), T the total exposure time, (^7) the 
differential size distribution (in terms of surface density), v 



the typical sky-projection relative speed, and d the typical 
distance to the TNOs. 

To derive (757) from the event rate, a functional form of 
the distribution needs to be assumed. On the other hand, if 
the integration is only over a small range of the size, we may 
derive an average value at that size. Noting that ^ 
^ s In 10, we have 



log s 



sv ds 



and 



dn 
dlog s 



dn 
d log s 



d 2 



31<S<S2 



«(S2 — Si) 

hill) 



V{S2 - S1 ) T lnl °W 



(9) 



(10) 



Assuming a typical distance d — 43 AU and a typical 
relative sky- project ion speed v — 30 km/s, with T — 240 ks 
and setting one detection as the upper limit, we have in the 
size range from 60 m to 300 m 



-P- ) < 1.7 x 10 10 deg" 
dlogs ' 



(11) 



The ecliptic latitude distribution of TNOs is not yet 
well determined. To be more precise, one should discuss 
the latitude with respect to the Kuiper Belt Plane (KBP) 
jElliot et al.ll2005l ) rather than to the ecliptic, but the differ- 
ence is not large. Sco X-l is 5.5° nort h of the ecliptic and 4.9° 
noth of the KBP. The CFHT survey (jTruiillo. Jewitt fc Luul 
l200ll) reported a 20° half-angle of an assumed gaussian dis- 
tribution in orbital inclination (with respect to the ecliptic). 
This half-angle in the inclination distribution translates to 
about 12.8° for a corresponding half-angle in the apparent 
ecliptic latitude distribution, assuming circular orbits. With 
this distribution, a factor of 1.14 should be applied to con- 
vert the estimate at the Sco X-l latitude to the e cliptic. How- 
ever, if we adopt the distribution reported in lElliot et al.l 
(2005), a factor of 4 should be applied for conversion to 
the KBP. In this pape r we u se the latitude distribution of 
iTruiillo. Jewitt fc Luul |200ll ). 



5.2 TNOs larger than 100 km 

For large TNOs, we use a rec ent result reported in 
iFuentes. George fc Holman] (|2009h . in which a double- 
power-law model describes well the luminosity function 
down to i?-band magnitude of about 27. With the notation 
employed in that paper, the differential luminosity function 
reads 



a(R) = C(10 
and 



-ojiffl-23) 



+10 



(a 2 - ai )(i? eq -23) 



10" 



2(H-23h 



" X (12) 



C = 0-23(1 + l0 (a2 -" l)(fl ^- 23) ) , (13) 

where a(R) is, with our notation here, equal to — The 
reported best fit values of the parameters in that distribution 



A rrO + 0.08 

are Qi = 0.73I 09 , a 2 



20+ 012 
u.^u_ 14 , 



i? cq = 25.01^, and 



023 



1.46!q'^2- C ' s therefore equal to 1.59. To convert the 



luminosity function into the size distribution, we assume a 
constant albedo and the same distance for all TNOs, and 
note that 



10 



-(- 



) . 



(14) 
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It follows that 
dn 
ds" 



C; 1 Ft— R e 



In 10 s c 



5 C 
In 10 s 



10 



ai(-Req-23) 



-1 
(15) 



and 
dn 
d log s 



dn 



s In 10 



ds 

5C , 1()ai (fl oq -23) 



+ 



(16) 



Adopting s cq = 90 km at R aq — 25, which as- 
sumes an albedo of 4% and a distance of 42 AU 
and adopts m_R = —27.6 for the J? - band magnitude 
of the Sun jFuentes. George fc Holmanl (|2009l ) . see also 
iFraser fc Kavelaarsl ^OOgD ), we have, at s = 90 km, that 



dn = 1.1 x 10 2 deg" 2 



d log s 



(17) 



5.3 TNOs of size at 500 m (observations of 
HST/FGS) 

ISchlichting et all {2009) reported the detection of an occul- 
tation event by a TNO of 500 m diameter in the HST/FGS 
archival data. They suggest a power index q = 3.9 ± 0.03 
for the TNO differential size distribution below 90 km and 
the accumulated number of TNOs larger than 500 m is 
2. ill * x 10 7 deg -2 . These values can be converted into 



(d 



log J 



in the following way. The accumulated distribution 



Ci 



■so 



1-8 



(18) 



and 
dn 
dlogs 



= ^slnlO 
ds 

= (l-q)C\ ( — Y 9 — In 10 



(19) 



For s = 500 m, C\ = 
have at s = 500 m that 

dn 



2.1±\ B 7 x 10 7 deg -2 . We therefore 



d log s 



= 1.4+i 2 x 10 8 deg" 



(20) 




0.010 0.100 1.000 10.000 100.000 1000. 000 

diameter (km) 



Figure 5. Size distributions of main-belt asteroids (MBAs) and 
trans-Neptunian objects (TNOs). Plotted here is the differential 
sky surface density at the ecliptic per decade of size in units 
of number per square degree. For TNOs, the upper limit at 
about 0.1 km, denoted with a downward arrow, is derived from 
our non-detection in the 240-ks RXTE/PCA data of Sco X-l, 
and is set at the level of one detection in 240 ks. The aster- 
isk symbol at 0.5 km is base d on the reported HST /FGS de- 
tection of an occultation event! Schlichting ct al. 2009). The thin 
solid curve is the double-power-la w distribution of large TNOs 
jFuentes. George fc Holmanl liooih . which shows a break at 90 
km. The HST/FGS result supports the existence of such a break, 
so does our current result. The thin dashed line is a direct extrap- 
olation from the large size end of the double-power-law towards 
smaller size and the thin dotted line is a power law anchoring on 
the double-power-law at 90 km with a power index of -3.0. We 
note that all the power indices in this figure, explicitly printed 
close to the corresponding power-law lines, are for the differential 
size distribution per decade of size, and are therefore the same as 
that of the cumulative size distribution. Upper limits derived from 
other ob servations are all above the thin dashed line and can be 
found in lSchlichting et alj j2009l V The double-power-law and the 
HST/FGS result in fact indicate a wavy shape for the TNO size 
distribution, which is similar to the behaviour of MBAs, plotted 
in the lower left half in this figure. The thick lines with power in- 
dices — 1.3 and —3.0 are for MBAs larger than 0.4 km, which arc 
directly observable. The MBA size distribution has a wavy shape 
with two bumps at about 100 km (not shown in this figure) and 5 
km. The thick dotted line with a power index of —2.5 is for a size 
distribution in collisional equilibrium assuming a constant struc- 
ture strength of MBAs in that size range. The thick dashed line 
with a power index of —5.8 is the one implied by Event El being 
caused by a 40-m MBA. The thick solid line with a power index 
of —4.0 is the upper limit at the level of assuming one detection 
down to 10 meters. 



6 SUMMARY AND DISCUSSION 

In the 240-ks RXTE/PCA data of Sco X-l taken from June 
2007 to October 2009, we found that millisecond dip events 
in the light curve are related to VLEs, which in turn suggests 
an instrumental origin for those millisecond dips. The rela- 
tion is strongest for Type-A VLEs and weakest for Type D. 
One dip event (denoted as Event El in this paper) at random 
probability of less than 7.5 x 10 -5 is found to be unrelated to 
any VLEs. The slight asymmetry, which is similar to many 
of the VLE-related dips, and the lacking of side lobes in the 
light curve of Event El cast some doubts on its nature be- 
ing an occultation event. However, since the count number 



is small (the central 4 bins in the light curve plotted in Fig- 
ure [2] have 2, 0, 5, and 2 counts respectively), one probably 
should treat the light curve with high uncertainty. It is also 
because of the low count number, diffraction pattern fitting 
does not provide strong enough information to constrain the 
property of the occulting body if it is indeed an occultation 
event, although acceptable fittings can be easily obtained. 
As discussed in previous sections, it might be due to a TNO 
of 150-m size, but with a rare retrograde orbit, or, it might 
be due to an MBA of 40-m size, but the associated detec- 
tion rate is incredibly high, or, it might be due to a very 
nearby object of meter size moving at a relative speed of a 
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few kilometers per second. Given all these uncertainties, we 
proceed to estimate the upper limits to the size distribution 
of TNOs and MBAs at the level of assuming one detection. 

The TNO size distribution has been stud- 
ied in several previous work s 

feoosh 



ISchlichting et all 

iFuentes. George fc Holmanf 



e.g.. | Bianco et al 



Bickerton, Kavclaar s fc Welch! 
Rogues et al l (|2006l ); 



(l2008h: iLiu et al l 



Fraser fc K avclaars 
20091) 



2008 



Bernstein et al.l l|2004h . In this paper 
we provide the most updated upper limit at the hectometer 
size, which is a size range that optical occultation surveys 
are unable to explore. Readers are referred to the caption 
of Figure [S] for more discussion. 

The MB A size distribution has been studied by sev- 
eral surveys l|lvezic et al. I l200ll ; iTedesco fc Desert! l2002i ; 
lYoshida fc Makamural 20071 ) down to the size of 400 me- 
ters. For smaller ones, if the collisional equilibrium and a 
constant tensile strength are assumed, a power law with 
power index -2.5 (cumulativ e) is expected (|Dohnanvill 19691 ; 
lO'Brien fc Greenberd I2OO3I ) . However, the tensile strength 
is known to be size dependent, and the size distribu- 
tion for the size range in the so-called strength-scaled 
regime in collisional equilibrium is expected to be steeper 
than the constant-strength solution (|Benz fc Asphauell 19991 ; 
lO'Brien fc Greenberd |2003T I . The upper limit that we esti- 
mate, as shown in Figure [5] has a power index —4.0, which 
in turn suggest s a power index of —2 .25 for the strength-size 
relation (jO'Brien fc Greenberd[2003r i . expressed in the form 
of the critical specific energy (conventionally denoted with 
Q* D ) as a function of size. This is still much steeper than all 
the reported strength-size relation in the literature to our 
knowledge, of which the steepest one has a p ower index of 
about -1 (|Durda. Greenberg fe Jedickdl 19981 ). 

Studying population properties of small TNOs and 
MBAs with serendipitous occultation events in X-rays is still 
in its infant stage, currently only possible for RXTE/PCA, 
the X-ray instrument with the largest effective area, ob- 
serving Sco X-l, the brightest X-ray source in the sky. It, 
however, suffers significantly from the instrumental dead 
time caused by high energy events. Observations with 
RXTE a nd future facilities such as ASTROSAT/LAXPC, 
AXTAR jChakrabartv et al.ll200Sl ), and IXO, together with 
a more rigorous design of t he survey, similar to those con- 
sider e d in the optical band (iBickerton. Welch fc Kavelaars I 
120091 ; iNihei et al.ll2007l ; lGaudill2004l ). will be able to yield 
more reliable and complete information, particularly when 
more X-ray sources in different directions can be employed. 
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